Evidence of hybridization between genetically distinct Baltic cod stocks during peak population abundance(s).

Range expansions can lead to increased contact of divergent populations, thus increasing the potential of hybridization events. Whether viable hybrids are produced will most likely depend on the level of genomic divergence and associated genomic incompatibilities between the different entities as well as environmental conditions. By taking advantage of historical Baltic cod (Gadus morhua) otolith samples combined with genotyping and whole genome sequencing, we here investigate the genetic impact of the increased spawning stock biomass of the eastern Baltic cod stock in the mid 1980s. The eastern Baltic cod is genetically highly differentiated from the adjacent western Baltic cod and locally adapted to the brackish environmental conditions in the deeper Eastern basins of the Baltic Sea unsuitable for its marine counterparts. Our genotyping results show an increased proportion of eastern Baltic cod in western Baltic areas (Mecklenburg Bay and Arkona Basin)-indicative of a range expansion westwards-during the peak population abundance in the 1980s. Additionally, we detect high frequencies of potential hybrids (including F1, F2 and backcrosses), verified by whole genome sequencing data for a subset of individuals. Analysis of mitochondrial genomes further indicates directional gene flow from eastern Baltic cod males to western Baltic cod females. Our findings unravel that increased overlap in distribution can promote hybridization between highly divergent populations and that the hybrids can be viable and survive under specific and favourable environmental conditions. However, the observed hybridization had seemingly no long-lasting impact on the continuous separation and genetic differentiation between the unique Baltic cod stocks.

Stable isotope analysis in food web research: Systematic review and a vision for the future for the Baltic Sea macro-region.

Food web research provides essential insights into ecosystem functioning, but practical applications in ecosystem-based management are hampered by a current lack of knowledge synthesis. To address this gap, we provide the first systematic review of ecological studies applying stable isotope analysis, a pivotal method in food web research, in the heavily anthropogenically impacted Baltic Sea macro-region. We identified a thriving research field, with 164 publications advancing a broad range of fundamental and applied research topics, but also found structural shortcomings limiting ecosystem-level understanding. We argue that enhanced collaboration and integration, including the systematic submission of Baltic Sea primary datasets to stable isotope databases, would help to overcome many of the current shortcomings, unify the scattered knowledge base, and promote future food web research and science-based resource management. The effort undertaken here demonstrates the value of macro-regional synthesis, in enhancing access to existing data and supporting strategic planning of research agendas.

Living apart together: Long-term coexistence of Baltic cod stocks associated with depth-specific habitat use.

Coexistence of fish populations (= stocks) of the same species is a common phenomenon. In the Baltic Sea, two genetically divergent stocks of Atlantic cod (Gadus morhua), Western Baltic cod (WBC) and Eastern Baltic cod (EBC), coexist in the Arkona Sea. Although the relative proportions of WBC and EBC in this area are considered in the current stock assessments, the mixing dynamics and ecological mechanisms underlying coexistence are not well understood. In this study, a genetically validated otolith shape analysis was used to develop the most comprehensive time series of annual stock mixing data (1977-2019) for WBC and EBC. Spatio-temporal mixing analysis confirmed that the two stocks coexist in the Arkona Sea, albeit with fluctuating mixing proportions over the 43-year observation period. Depth-stratified analysis revealed a strong correlation between capture depth and stock mixing patterns, with high proportions of WBC in shallower waters (48-61% in <20m) and increasing proportions of EBC in deeper waters (50-86% in 40-70m). Consistent depth-specific mixing patterns indicate stable differences in depth distribution and habitat use of WBC and EBC that may thus underlie the long-term coexistence of the two stocks in the Arkona Sea. These differences were also reflected in significantly different proportions of WBC and EBC in fisheries applying passive gears in shallower waters (more WBC) and active gears in deeper waters (more EBC). This highlights the potential for fishing gear-specific exploitation of different stocks, and calls for stronger consideration of capture depth and gear type in stock assessments. This novel evidence provides the basis for improved approaches to research, monitoring and management of Baltic cod stocks.

Niche partitioning between planktivorous fish in the pelagic Baltic Sea assessed by DNA metabarcoding, qPCR and microscopy.

Marine communities undergo rapid changes related to human-induced ecosystem pressures. The Baltic Sea pelagic food web has experienced several regime shifts during the past century, resulting in a system where competition between the dominant planktivorous mesopredatory clupeid fish species herring (Clupea harengus) and sprat (Sprattus sprattus) and the rapidly increasing stickleback (Gasterosteus aculeatus) population is assumed to be high. Here, we investigate diet overlap between these three planktivorous fishes in the Baltic Sea, utilizing DNA metabarcoding on the 18S rRNA gene and the COI gene, targeted qPCR, and microscopy. Our results show niche differentiation between clupeids and stickleback, and highlight that rotifers play an important role in this pattern, as a resource that is not being used by the clupeids nor by other zooplankton in spring. We further show that all the diet assessment methods used in this study are consistent, but also that DNA metabarcoding describes the plankton-fish link at the highest taxonomic resolution. This study suggests that rotifers and other understudied soft-bodied prey may have an important function in the pelagic food web and that the growing population of pelagic stickleback may be supported by the open feeding niche offered by the rotifers.

Food web assessments in the Baltic Sea: Models bridging the gap between indicators and policy needs.

Ecosystem-based management requires understanding of food webs. Consequently, assessment of food web status is mandatory according to the European Union's Marine Strategy Framework Directive (MSFD) for EU Member States. However, how to best monitor and assess food webs in practise has proven a challenging question. Here, we review and assess the current status of food web indicators and food web models, and discuss whether the models can help addressing current shortcomings of indicator-based food web assessments, using the Baltic Sea as an example region. We show that although the MSFD food web assessment was designed to use food web indicators alone, they are currently poorly fit for the purpose, because they lack interconnectivity of trophic guilds. We then argue that the multiple food web models published for this region have a high potential to provide additional coherence to the definition of good environmental status, the evaluation of uncertainties, and estimates for unsampled indicator values, but we also identify current limitations that stand in the way of more formal implementation of this approach. We close with a discussion of which current models have the best capacity for this purpose in the Baltic Sea, and of the way forward towards the combination of measurable indicators and modelling approaches in food web assessments.

DNA metabarcoding reveals the importance of gelatinous zooplankton in the diet of Pandalus borealis, a keystone species in the Arctic.

Information about the dietary composition of a species is crucial to understanding their position and role in the food web. Increasingly, molecular approaches such as DNA metabarcoding are used in studying trophic relationships, not least because they may alleviate problems such as low taxonomic resolution or underestimation of digestible taxa in the diet. Here, we used DNA metabarcoding with universal primers for cytochrome c oxidase I (COI) to study the diet composition of the northern shrimp (Pandalus borealis), an Arctic keystone species with large socio-economic importance. Across locations, jellyfish and chaetognaths were the most important components in the diet of P. borealis, jointly accounting for 40%-60% of the total read abundance. This dietary importance of gelatinous zooplankton contrasts sharply with published results based on stomach content analysis. At the same time, diet composition differed between fjord and shelf locations, pointing to different food webs supporting P. borealis in these two systems. Our study underlines the potential of molecular approaches to provide new insights into the diet of marine invertebrates that are difficult to obtain with traditional methods, and calls for a revision of the role of gelatinous zooplankton in the diet of the key Arctic species P. borealis, and in extension, Arctic food webs.

Multi-trophic markers illuminate the understanding of the functioning of a remote, low coral cover Marquesan coral reef food web.

We studied the food web structure and functioning of a coral reef ecosystem in the Marquesas Islands, French Polynesia, characterized by low coral cover, high sea surface temperature and meso- to eutrophic waters. The Marquesas constitute a relevant ecosystem to understand the functioning of low diversity reefs that are also subject to global change. A multi-tracer assessment of organic matter pathways was run to delineate ecosystem functioning, using analysis of fatty acids, bulk and compound specific stable isotope analysis and stable isotopes mixing models. Macroalgae and phytoplankton were the two major food sources fueling this food web with, however, some marked seasonal variations. Specifically, zooplankton relied on phytoplankton-derived organic matter and herbivorous fishes on macroalgae-derived organic matter to a much higher extent in summer than in winter (~ 75% vs. ~ 15%, and ~ 70 to 75% vs. ~ 5 to 15%, respectively) . Despite remarkably high δ15N values for all trophic compartments, likely due to local dynamics in the nitrogen stock, trophic levels of consumers were similar to those of other coral reef ecosystems. These findings shed light on the functioning of low coral cover systems, which are expected to expand worldwide under global change.

Characterizing niche differentiation among marine consumers with amino acid δ13C fingerprinting.

Marine food webs are highly compartmentalized, and characterizing the trophic niches among consumers is important for predicting how impact from human activities affects the structuring and functioning of marine food webs. Biomarkers such as bulk stable isotopes have proven to be powerful tools to elucidate trophic niches, but they may lack in resolution, particularly when spatiotemporal variability in a system is high. To close this gap, we investigated whether carbon isotope (δ13C) patterns of essential amino acids (EAAs), also termed δ13CAA fingerprints, can characterize niche differentiation in a highly dynamic marine system. Specifically, we tested the ability of δ13CAA fingerprints to differentiate trophic niches among six functional groups and ten individual species in the Baltic Sea. We also tested whether fingerprints of the common zooplanktivorous fishes, herring and sprat, differ among four Baltic Sea regions with different biochemical conditions and phytoplankton assemblages. Additionally, we investigated how these results compared to bulk C and N isotope data for the same sample set. We found significantly different δ13CAA fingerprints among all six functional groups. Species differentiation was in comparison less distinct, due to partial convergence of the species' fingerprints within functional groups. Herring and sprat displayed region-specific δ13CAA fingerprints indicating that this approach could be used as a migratory marker. Niche metrics analyses showed that bulk isotope data had a lower power to differentiate between trophic niches than δ13CAA fingerprinting. We conclude that δ13CAA fingerprinting has a strong potential to advance our understanding of ecological niches, and trophic linkages from producers to higher trophic levels in dynamic marine systems. Given how management practices of marine resources and habitats are reshaping the structure and function of marine food webs, implementing new and powerful tracer methods are urgently needed to improve the knowledge base for policy makers.

Between source and sea: The role of wastewater treatment in reducing marine microplastics.

Wastewater treatment plants (WWTPs) are a focal point for the removal of microplastic (MP) particles before they are discharged into aquatic environments. WWTPs are capable of removing substantial quantities of larger MP particles but are inefficient in removing particles with any one dimension of less than 100 µm, with influents and effluents tending to have similar quantities of these smaller particles. As a single WWTP may release >100 billion MP particles annually, collectively WWTPs are significant contributors to the problem of MP pollution of global surface waters. Currently, there are no policies or regulations requiring the removal of MPs during wastewater treatment, but as concern about MP pollution grows, the potential for wastewater technologies to capture particles before they reach surface waters has begun to attract attention. There are promising technologies in various stages of development that may improve the removal of MP particles from wastewater. Better incentivization could speed up the research, development and adoption of innovative practices. This paper describes the current state of knowledge regarding MPs, wastewater and relevant policies that could influence the development and deployment of new technologies within WWTPs. We review existing technologies for capturing very small MP particles and examine new developments that may have the potential to overcome the shortcomings of existing methods. The types of collaborations needed to encourage and incentivize innovation within the wastewater sector are also discussed, specifically strong partnerships among scientific and engineering researchers, industry stakeholders, and policy decision makers.

Assessing SNP-markers to study population mixing and ecological adaptation in Baltic cod.

Atlantic cod (Gadus morhua) is a species of great ecological and economical importance in the Baltic Sea. Here, two genetically differentiated stocks, the western and the eastern Baltic cod, display substantial mechanical mixing, hampering our understanding of cod ecology and impeding stock assessments and management. Based on whole-genome re-sequencing data from reference samples obtained from the study area, we designed two different panels of Single Nucleotide Polymorphisms markers (SNPs), which take into account the exceptional genome architecture of cod. A minimum panel of 20 diagnostic SNPs and an extended panel (20 diagnostic and 18 biologically informative SNPs, 38 in total) were developed and validated to distinguish unambiguously between the western and the eastern Baltic cod stocks and to enable studies of local adaptation to the specific environment in the Baltic Sea, respectively. We tested both panels on cod sampled from the southern Baltic Sea (n = 603) caught in 2015 and 2016. Genotyping results showed that catches from the mixing zone in the Arkona Sea, were composed of similar proportions of individuals of the western and the eastern stock. Catches from adjacent areas to the east, the Bornholm Basin and Gdansk Deep, were exclusively composed of eastern Baltic cod, whereas catches from adjacent western areas (Belt Sea and Öresund) were composed of western Baltic cod. Interestingly, the two Baltic cod stocks showed strong genetic differences at loci associated with life-history trait candidate genes, highlighting the species' potential for ecological adaptation even at small geographical scales. The minimum and the extended panel of SNP markers presented in this study provide powerful tools for future applications in research and fisheries management to further illuminate the mixing dynamics of cod in the Baltic Sea and to better understand Baltic cod ecology.

Genetic analyses reveal complex dynamics within a marine fish management area.

Genetic data have great potential for improving fisheries management by identifying the fundamental management units-that is, the biological populations-and their mixing. However, so far, the number of practical cases of marine fisheries management using genetics has been limited. Here, we used Atlantic cod in the Baltic Sea to demonstrate the applicability of genetics to a complex management scenario involving mixing of two genetically divergent populations. Specifically, we addressed several assumptions used in the current assessment of the two populations. Through analysis of 483 single nucleotide polymorphisms (SNPs) distributed across the Atlantic cod genome, we confirmed that a model of mechanical mixing, rather than hybridization and introgression, best explained the pattern of genetic differentiation. Thus, the fishery is best monitored as a mixed-stock fishery. Next, we developed a targeted panel of 39 SNPs with high statistical power for identifying population of origin and analyzed more than 2,000 tissue samples collected between 2011 and 2015 as well as 260 otoliths collected in 2003/2004. These data provided high spatial resolution and allowed us to investigate geographical trends in mixing, to compare patterns for different life stages and to investigate temporal trends in mixing. We found similar geographical trends for the two time points represented by tissue and otolith samples and that a recently implemented geographical management separation of the two populations provided a relatively close match to their distributions. In contrast to the current assumption, we found that patterns of mixing differed between juveniles and adults, a signal likely linked to the different reproductive dynamics of the two populations. Collectively, our data confirm that genetics is an operational tool for complex fisheries management applications. We recommend focussing on developing population assessment models and fisheries management frameworks to capitalize fully on the additional information offered by genetically assisted fisheries monitoring.

Disentangling structural genomic and behavioural barriers in a sea of connectivity.

Genetic divergence among populations arises through natural selection or drift and is counteracted by connectivity and gene flow. In sympatric populations, isolating mechanisms are thus needed to limit the homogenizing effects of gene flow to allow for adaptation and speciation. Chromosomal inversions act as an important mechanism maintaining isolating barriers, yet their role in sympatric populations and divergence with gene flow is not entirely understood. Here, we revisit the question of whether inversions play a role in the divergence of connected populations of the marine fish Atlantic cod (Gadus morhua), by exploring a unique data set combining whole-genome sequencing data and behavioural data obtained with acoustic telemetry. Within a confined fjord environment, we find three genetically differentiated Atlantic cod types belonging to the oceanic North Sea population, the western Baltic population and a local fjord-type cod. Continuous behavioural tracking over 4 year revealed temporally stable sympatry of these types within the fjord. Despite overall weak genetic differentiation consistent with high levels of gene flow, we detected significant frequency shifts of three previously identified inversions, indicating an adaptive barrier to gene flow. In addition, behavioural data indicated that North Sea cod and individuals homozygous for the LG12 inversion had lower fitness in the fjord environment. However, North Sea and fjord-type cod also occupy different depths, possibly contributing to prezygotic reproductive isolation and representing a behavioural barrier to gene flow. Our results provide the first insights into a complex interplay of genomic and behavioural isolating barriers in Atlantic cod and establish a new model system towards an understanding of the role of genomic structural variants in adaptation and diversification.

A fundamental difference between macrobiota and microbial eukaryotes: protistan plankton has a species maximum in the freshwater-marine transition zone of the Baltic Sea.

Remane's Artenminimum at the horohalinicum is a fundamental concept in ecology to describe and explain the distribution of organisms along salinity gradients. However, a recent metadata analysis challenged this concept for protists, proposing a species maximum in brackish waters. Due to data bias, this literature-based investigation was highly discussed. Reliable data verifying or rejecting the species minimum for protists in brackish waters were critically lacking. Here, we sampled a pronounced salinity gradient along a west-east transect in the Baltic Sea and analysed protistan plankton communities using high-throughput eDNA metabarcoding. A strong salinity barrier at the upper limit of the horohalinicum and 10 psu appeared to select for significant shifts in protistan community structures, with dinoflagellates being dominant at lower salinities, and dictyochophytes and diatoms being keyplayers at higher salinities. Also in vertical water column gradients in deeper basins (Kiel Bight, Arkona and Bornholm Basin) appeared salinity as significant environmental determinant influencing alpha- and beta-diversity patterns. Importantly, alpha-diversity indices revealed species maxima in brackish waters, that is, indeed contrasting Remane's Artenminimum concept. Statistical analyses confirmed salinity as the major driving force for protistan community structuring with high significance. This suggests that macrobiota and microbial eukaryotes follow fundamentally different rules regarding diversity patterns in the transition zone from freshwater to marine waters.

The Baltic Sea as a time machine for the future coastal ocean.

Coastal global oceans are expected to undergo drastic changes driven by climate change and increasing anthropogenic pressures in coming decades. Predicting specific future conditions and assessing the best management strategies to maintain ecosystem integrity and sustainable resource use are difficult, because of multiple interacting pressures, uncertain projections, and a lack of test cases for management. We argue that the Baltic Sea can serve as a time machine to study consequences and mitigation of future coastal perturbations, due to its unique combination of an early history of multistressor disturbance and ecosystem deterioration and early implementation of cross-border environmental management to address these problems. The Baltic Sea also stands out in providing a strong scientific foundation and accessibility to long-term data series that provide a unique opportunity to assess the efficacy of management actions to address the breakdown of ecosystem functions. Trend reversals such as the return of top predators, recovering fish stocks, and reduced input of nutrient and harmful substances could be achieved only by implementing an international, cooperative governance structure transcending its complex multistate policy setting, with integrated management of watershed and sea. The Baltic Sea also demonstrates how rapidly progressing global pressures, particularly warming of Baltic waters and the surrounding catchment area, can offset the efficacy of current management approaches. This situation calls for management that is (i) conservative to provide a buffer against regionally unmanageable global perturbations, (ii) adaptive to react to new management challenges, and, ultimately, (iii) multisectorial and integrative to address conflicts associated with economic trade-offs.

No increase in marine microplastic concentration over the last three decades - A case study from the Baltic Sea.

Microplastic is considered a potential threat to marine life as it is ingested by a wide variety of species. Most studies on microplastic ingestion are short-term investigations and little is currently known about how this potential threat has developed over the last decades where global plastic production has increased exponentially. Here we present the first long-term study on microplastic in the marine environment, covering three decades from 1987 to 2015, based on a unique sample set originally collected and conserved for food web studies. We investigated the microplastic concentration in plankton samples and in digestive tracts of two economically and ecologically important planktivorous forage fish species, Atlantic herring (Clupea harengus) and European sprat (Sprattus sprattus), in the Baltic Sea, an ecosystem which is under high anthropogenic pressure and has undergone considerable changes over the past decades. Surprisingly, neither the concentration of microplastic in the plankton samples nor in the digestive tracts changed significantly over the investigated time period. Average microplastic concentration in the plankton samples was 0.21±0.15particlesm-3. Of 814 fish examined, 20% contained plastic particles, of which 95% were characterized as microplastic (<5mm) and of these 93% were fibres. There were no significant differences in the plastic content between species, locations, or time of day the fish were caught. However, fish size and microplastic in the digestive tracts were positively correlated, and the fish contained more plastic during summer than during spring, which may be explained by increased food uptake with size and seasonal differences in feeding activity. This study highlights that even though microplastic has been present in the Baltic environment and the digestive tracts of fishes for decades, the levels have not changed in this period. This underscores the need for greater understanding of how plastic is cycled through marine ecosystems. The stability of plastic concentration and contamination over time observed here indicates that the type and level of microplastic pollution may be more closely correlated to specific human activities in a region than to global plastic production and utilization as such.

Ancient DNA reveals the Arctic origin of Viking Age cod from Haithabu, Germany.

Knowledge of the range and chronology of historic trade and long-distance transport of natural resources is essential for determining the impacts of past human activities on marine environments. However, the specific biological sources of imported fauna are often difficult to identify, in particular if species have a wide spatial distribution and lack clear osteological or isotopic differentiation between populations. Here, we report that ancient fish-bone remains, despite being porous, brittle, and light, provide an excellent source of endogenous DNA (15-46%) of sufficient quality for whole-genome reconstruction. By comparing ancient sequence data to that of modern specimens, we determine the biological origin of 15 Viking Age (800-1066 CE) and subsequent medieval (1066-1280 CE) Atlantic cod (Gadus morhua) specimens from excavation sites in Germany, Norway, and the United Kingdom. Archaeological context indicates that one of these sites was a fishing settlement for the procurement of local catches, whereas the other localities were centers of trade. Fish from the trade sites show a mixed ancestry and are statistically differentiated from local fish populations. Moreover, Viking Age samples from Haithabu, Germany, are traced back to the North East Arctic Atlantic cod population that has supported the Lofoten fisheries of Norway for centuries. Our results resolve a long-standing controversial hypothesis and indicate that the marine resources of the North Atlantic Ocean were used to sustain an international demand for protein as far back as the Viking Age.

Temporal dietary shift in jellyfish revealed by stable isotope analysis.

A temporal change in the stable isotope (SI) composition of jellyfish in the Kiel Fjord, Western Baltic Sea, was documented by analyzing δ13C, δ15N and δ34S of bell tissue of Aurelia aurita and Cyanea capillata in the period between June and October 2011. A strong and significant temporal change in all SI values of A. aurita was found, including an increase of ~3 ‰ in δ13C, a decrease of ~4 ‰ in δ15N and sharp decline of ~7 ‰ in δ34S. While knowledge gaps in jellyfish isotope ecology, in particular the lack of reliable trophic enrichment factors, call for a conservative interpretation of our data, observed changes in particular in δ34S, as indicated by means of a MixSIR mixing model, would be consistent with a temporal dietary shift in A. aurita from mesozooplankton (>150 µm) to microplankton and small re-suspended particles (0.8-20 µm) from the benthos. Presence of a hitherto unidentified food source not included in the model could also contribute to the shift. During the 2-month occurrence of C. capillata, its isotope composition remained stable and was consistent with a mainly mesozooplanktonic diet. Mixing model output, mainly driven by δ34S values, indicated a lower proportion of A. aurita in the diet of C. capillata than previously reported, and thus to a potentially lesser importance of intraguild predation among jellyfish in the Kiel Fjord. Overall, our results clearly highlighted the potential for substantial intraspecific isotopic seasonal variation in jellyfish, which should be taken into account in future feeding ecology studies on this group.

The great melting pot. Common sole population connectivity assessed by otolith and water fingerprints.

Quantifying the scale and importance of individual dispersion between populations and life stages is a key challenge in marine ecology. The common sole (Solea solea), an important commercial flatfish in the North Sea, Atlantic Ocean and the Mediterranean Sea, has a marine pelagic larval stage, a benthic juvenile stage in coastal nurseries (lagoons, estuaries or shallow marine areas) and a benthic adult stage in deeper marine waters on the continental shelf. To date, the ecological connectivity among these life stages has been little assessed in the Mediterranean. Here, such an assessment is provided for the first time for the Gulf of Lions, NW Mediterranean, based on a dataset on otolith microchemistry and stable isotopic composition as indicators of the water masses inhabited by individual fish. Specifically, otolith Ba/Ca and Sr/Ca profiles, and δ(13)C and δ(18)O values of adults collected in four areas of the Gulf of Lions were compared with those of young-of-the-year collected in different coastal nurseries. Results showed that a high proportion of adults (>46%) were influenced by river inputs during their larval stage. Furthermore Sr/Ca ratios and the otolith length at one year of age revealed that most adults (∼70%) spent their juvenile stage in nurseries with high salinity, whereas the remainder used brackish environments. In total, data were consistent with the use of six nursery types, three with high salinity (marine areas and two types of highly saline lagoons) and three brackish (coastal areas near river mouths, and two types of brackish environments), all of which contributed to the replenishment of adult populations. These finding implicated panmixia in sole population in the Gulf of Lions and claimed for a habitat integrated management of fisheries.

Anthropogenic hybridization between endangered migratory and commercially harvested stationary whitefish taxa (Coregonus spp.).

Natural hybridization plays a key role in the process of speciation. However, anthropogenic (human induced) hybridization of historically isolated taxa raises conservation issues. Due to weak barriers to gene flow and the presence of endangered taxa, the whitefish species complex is an excellent study system to investigate the consequences of hybridization in conservation. We focused on three naturally reproductively isolated whitefish taxa in Germany: the endangered, anadromous North Sea houting (NSH) and Baltic houting (BH), which were reintroduced after local extinction, and the commercially stocked European whitefish (EW). To evaluate the genetic integrity of each taxon, source and reintroduced populations of NSH and BH, and EW populations were characterized based on two mitochondrial and 17 microsatellite loci. Additionally, we investigated gill raker counts as an adaptive phenotypic trait. Even though clear genetic and phenotypic differentiation confirmed the houtings as separate evolutionarily significant units, admixture analyses revealed an extensive hybrid zone. Hybridizations were introgressive, positively correlated with genetic diversity, and were reflected in the gill raker counts. The BH distribution range showed higher heterogeneity and stronger admixture than the NSH range. Erroneous stocking with non-native genotypes best explained these patterns, which pose challenges for the conservation of the endangered NSH and BH.